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We explore the zero temperature phase behavior of a two-dimensional two-component atomic 
Fermi gas with population and mass imbalance in the regime of the BEC-BCS crossover. Working 
in the mean-field approximation, we show that the normal and homogeneous balanced superfluid 
phases are separated by an inhomogeneous superfluid phase of Fulde-Ferrel-Larkin-Ovchinnikov 
(FFLO) type. We obtain an analytical expression for the line of continuous transitions separating 
the normal and inhomogeneous FFLO phases. We further show that the transition from the FFLO 
phase to the homogeneous balanced superfluid is discontinuous leading to phase separation. If the 
species have different masses, the superfluid phase is favored when the lighter species is in excess. 
We explore the implications of these findings for the properties of the two-component Fermi gas in 
the atomic trap geometry. Finally, we compare and contrast our findings with the predicted phase 
behavior of the electron-hole bilayer system. 
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I. INTRODUCTION 



By controlling interaction through a magnetically- 
tuned Feshbach resonance, ultracold atomic Fermi gases 
have provided a versatile arena in which to explore pair- 
ing phenomena and superfluidity [J, [3, 3 ■ Already the 
crossover between the Bose-Einstein condensate (BEG) 
phase of strongly bound diatomic molecules to the BGS 
phase of weakly bound Gooper pairs has been observed 
experimentally [l|, H, H, [y, 0, §• In recent years, 
much attention has been focused on the phase behav- 
ior of two-component Fermi gases with pop ulation im- 

WM, MM, M M,\m M M HH, 

\2l\ . and generalized mass ratios be- 
tween different species [H, M, M,M,M,M,M,M,^- 
The symmetry breaking effect of population and mass 
imbalance destabilizes the condensate leading to an en- 
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riched phase diagram characterized by tricritical point 
behavior with first order transitions separating normal 
and superfluid phases at low temperatures [36|. More 
detailed studies have shown that, on the weak cou- 
pling side of the crossover, the transition into a homo- 
geneous superfluid phase at low temperatures is pre- 
empted by the dev elop ment of an inhomogeneous super- 
fluid phase d, M, M, [13, m, m, EO, m\M, Si- This 
is a manifestation of the FFLO phase predicted to occur 



in superconducting electron systems subject to a Zeeman 
fleld |44l. [dsl . In the three-dimensional system, the FFLO 
phase is predicted to occupy only a small region of the 
phase diagram making its experimental identification in 
the atomic trap geometry challenging. Indeed, even in 
solid state systems, the FFLO phase has only recently 
been observed [46j . 

The potential for an FFLO instability at a single wave 
vector in a three-dimensional ultracold atomic gas with 
only population imbalance was explored by Hu and Liu 
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jl7j and Zhang and Duan [2J|. They found a stable 
FFLO phase only on the BGS side of the resonance. Ad- 
ditionally, Wu and Yip [13] showed the three-dimensional 
system is unstable to FFLO superfiuid currents, but these 
were not found in the non-uniform three-dimensional 
trap experiments of Shin et al. [ist . 

Lately, efforts have been made to explore the effects 
of population imbalance on pairing in two-component 
Fermi gases in two-dimensions. Although the phase dia- 
gram of the zero temperature system has been explored 
in the regime of BEC-BGS crossover in the mean-field 
approximation [2^, the potential for FFLO phase for- 
mation has not yet been addressed. By contrast, mo- 
tivated by potential applications to strongly anisotropic 
layered systems, several theoretical studies have explored 
the potential for superconducting FFLO phase forma- 
tion in two-dimensional electron systems [49l . [sol . 511. 
A quasiclassical analysis by Gombescot and Mora 52| . 
involving a Ginzburg-Landau expansion of the free en- 
ergy in Fourier components of the superconducting or- 
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der parameter, found that the FFLO transition in two 
dimensions was continuous at low temperatures. In a 
separate study of condensation in electron-hole bilayers, 
Fieri et al. [50| argued that the FFLO phase can occupy 
a significant part of the two-dimensional phase diagram. 
Motivated by these investigations, and the potential im- 
pact on the atomic gas system, in the following we will 
investigate the potential for FFLO phase formation in 
the two-dimensional two-component Fermi gas address- 
ing both population imbalance and generalized mass ra- 
tios. 

In the context of ultracold atomic Fermi gases, a two- 
dimensional system can be experimentally realized by 
confining the gas with a one-dimensional optical lat- 
tice consisting of two counter-propagating laser beams 
[ssl . [S^ . [55I . [Sq . These impose a periodic potential, with 
antinodes spaced every half wavelength. The interwell 
barrier energy, and therefore the tunneling rate, depends 
on the laser intensity, which should be chosen to be much 
higher than the species chemical potential and the super- 
fiuid gap [53, nil- This inhibits transfer of atoms between 
layers and the Fermi gas separates into stacked quasi two- 
dimensional layers. The depth of the optical potential 
can always be varied independently of the external har- 
monic trapping potential and species chemical potentials 
so it should always be possible to reduce the tunneling 
rate sufficiently that the cold atom gas can be regarded 
as being two-dimensional gas. Within a layer, the short- 
ranged interaction of the atoms can be adjusted by ex- 
ploiting a Feshbach resonance. It has been suggested [5§| 
that due to the possible formation of dressed molecules 
a single band theory could incorrectly predict cloud size 
in the strong coupling limit. However, here we are inter- 
ested in the weak coupling limit and the behavior at the 
BEC-BCS crossover where we believe that a single band 
theory will encompass the essential behavior allowing us 
to capture the qualitative phase structure. 

Finally, further motivation for the investigation of su- 
perfluidity in the mass imbalanced system comes from 
studies of cxciton condensation in bilayer electron-hole 
systems. In recent years, attempts to realize a condensed 
exciton phase have focused on quantum well structures 
where electrons and holes are restricted to neighbor- 
ing two-dimensional layers [i^, [s^, [lH. The range of 
the Coulomb interaction between the particles can be 
shortened by introducing a screening layer. As with 
the two-component Fermi gas, the electron-hole system 
affords the possibility of tuning between a supcrfiuid 
of tightly-bound pairs (excitons) to a condensate phase 
of an electron-hole plasma. Moreover, while one can, 
in principle, engineer a balanced electron-hole popula- 
tion, the effective masses of the electron and hole quasi- 
particles in the semiconductor are typically quite differ- 
ent. In GaAs, the ratio of the hole to electron mass is ca. 
iTLh/mc = 4.3. Condensation phenomena in mass imbal- 
anced systems have also been explored in the context of 
quantum chromo dynamics, where the particles represent 
different species of quarks [B^I ■ 



The remainder of the paper is organized as follows: In 
Sec. lUwe begin by deriving an expression for the ther- 
modynamic potential in the mean-field approximation for 
the two-component Fermi gas allowing for the develop- 
ment of an inhomogeneous condensate phase. In Sec. IIIII 
we use this result to elucidate the zero-temperature phase 
diagram of the system for a uniform order parameter 
both at fixed chemical potential and fixed number den- 
sity. In Sec. IIVI we explore the tendency of the system 
to condense into an inhomogeneous superfluid phase. In 
particular we combine the results of a Ginzburg-Landau 
expansion with the numerical analysis of the thermody- 
namic potential to infer the region over which the inho- 
mogeneous phase persists. Finally, in Scc.|Vl we examine 
the properties of the atomic Fermi gas in the harmonic 
trap geometry, concluding our discussion in Sec. IVII 



II. MEAN-FIELD THEORY 

Our starting point is a two-component Fermi gas with 
each species indexed by a pseudo-spin a e {Tii} = 
{+1,-1}. The single-particle dispersion ek,o- = fc^/2mo- 
depends on the different effective masses of the two 
species rricr, throughout we set h = 1. Introducing the 
reduced mass, 1/mR = (l/frif -I- \/m{)/2 and the mass 
ratio r = m^/m-^ we have m-f = tor(1 -|- l/r)/2 and 
my — 771^(1 + r)/2. To enforce a population imbalance, 
each species must be associated with its own chemical 
potential, fi^ = ^ + ah. With these definitions, one may 
see that the phase diagram of the system is symmetric 
under the transformation, {h,r) t-^ {—h, l/r). 

In the following, we will focus on the zero temperature 
phase behavior of the system as predicted by mean-field 
theory. In doing so, we will miss non-perturbative effects 
that appear at large mass ratios. In particular, when the 
ratio of masses is greater than 13.6, it is known that, in 
three-dimensions, two heavy and one light fermion can 
form a three-body weakly bound state [6l|, [H, [g^ . Our 
analysis does not include the possibility of such a state. 

To explore the regime of BEC-BCS crossover, we will 
focus our attention on a single-channel Hamiltonian de- 
scribing a wide Feshbach resonance where the closed 
channel population remains small throughout [2^. [6^. [65| . 
The quantum partition function for the system can then 
be expressed as a functional field integral over fermionic 
fields -0^ and -0^, Z = ^ e'^^'^'^'^^'Dil^'DTp, with the action 



dr 



^"^'k.cr {dr +Ck,o-)'(/'k., 



-\ X! V'k,T^q-k,J,^k'-kV'k'aV'q-k', 



k,k'.q 



where V denotes the two-body interaction potential, and 
(3 = l/fceT is the inverse temperature. Here, for brevity, 
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we have set ^k,o- = ek.o- — Mcr- Anticipating the devel- channel, with Z = J e~^^'^'^'^'^^'DipVtp'DAVA, where 
opment of pair correlations, we introduce a Hubbard- the action now takes the form, 
Stratonovich decoupling of the interaction in the Cooper 



Lj,k,k' 

j-^^ ^ ^' ^' ' \ Ao,q -lUJ - 4k+q/2,i / V '/'<^,q/2+k,i / 

I 



Anticipating that the transition to the supcrfluid (SF) 
from the normal phase is continuous (a property al- 
ready established in the weak coupling limit of the two- 
dimensional system by Combcscot and Mora (s^). we 
will suppose that the order parameter is characterized by 
a single plane- wave state corresponding to the statio nary 
saddle-point solution, A^^^q ~ AQ(5q_Q(5tj,o [13, m, m- 
In this case, Q = describes the homogeneous SF state 
while, for Q ^ 0, the condensate is of FFLO type. If the 
transition to the inhomogeneous phase is first order, the 
single wave vector assumption necessitates some degree 
of approximation that will underestimate the width of 
the FFLO region in the phase diagram. 

Then, approximating the functional integral over fields 
A by its mean-field value, and taking the interaction to 
be contact, V{r) = g6^{r), integration over the fermionic 
fields gives the thermodynamic potential. 



^Ql 



t 



--Trln 



1 



(1) 



where Ck,± = (Ck-Q/2,T ± Ck+q/2,i)/2, and Ei, = (^g _^ 

■om this express: 
larization or "magnetization" 



AqP)^/^. From this expression, one can obtain the po- 



rn = n| — n| = —dQ/dh 

= riF [Ek - Ck,-) - riF (-Bk + Ck,-) , (2) 

and the total number density, 
n = n| + ?i| = — dr2/d/.j 

^""'^ {UF (-Bk - Ck,-) + np (^k + Ck,-) - 
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where nF{E) ~ l/(l-|-e~^'^) denotes the Fermi function. 

Finally, to regularize the unphysical UV divergences 
associated with the 5-function form of contact interaction 
(and contained within the term labeled by | in Eqn. ([1])), 
we will set 
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g Jo '2E + 



-dE. 



(4) 



where E^ denotes the energy of the two-body bound 
state m, 1131 ■ will then be used as a control param- 
eter to tune through the BEC-BCS crossover. As E'b is 
increased, the system evolves continuously from the weak 
coupling BCS phase to the strong coupling BEC phase 
of tightly-bound pairs. 

Having obtained the thermodynamic potential in the 
mean-field approximation, we now outline our strategy 
for calculating the zero temperature phase diagram. As 
a platform to address the potential for inhomogeneous 
phase formation, in the following section we begin by 
establishing the phase diagram associated with a uni- 
form order parameter, i.e. Q = 0. In this case, the 
integrations associated with the thermodynamic poten- 
tial Eqn. ([T]) can be evaluated analytically and many key 
features of the generalized phase diagram understood. 
Then, in Sec. lIVi we return to the more general situation, 
exploring the capacity for inhomogeneous phase forma- 
tion. After confirming that, in the single wave vector 
approximation, the transition to the SF phase is always 
continuous, we develop a Ginzburg-Landau expansion of 
the action to determine the locus of the normal-FFLO 
phase boundary analytically. Combining these results, 
we determine the phase diagram for a spatially uniform 
system as function of fixed chemical potential and, sep- 
arately, as a function of fixed particle number. Finally, 
in Sec. |Vl we apply these results to the problem of reso- 
nance superfluidity in the physically realizable harmoni- 
cally trapped system. 



III. UNIFORM SUPERFLUID 



In this section we focus on the phase diagram of a 
system in which the order parameter is constrained to 
be uniform. At zero temperature, setting Q = 0, the 
thermodynamic potential can be evaluated analytically 
for arbitrary population imbalance and mass ratio. 
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-0 /i'^-IAol 



X \^\h'\{fi, + 2fif) - fi^^i'^ + |AoP - |Ao|' In 
+ sgn{h') (^Ri^,^-h'^) + \^Hr)^^ 



(5) 



Here denotes the Heavisidc step function, = TOR/27r 
the two-dimensional density of states of the reduced mass 
system and, for clarity, we have defined the set of param- 
eters, 



R: 



^J■c 



r-l 



h' = 



h — iiR 



Vh'^-\Ao\ 



= AtVl - ■ (6) 



By minimizing the thermodynamic potential with respect 
to Ao, one obtains the loci of phase boundaries shown in 
Table m When the mass ratio is unity {R = 0), these 
results coincide with those obtained in Ref . [1^ . In par- 
ticular, one may note that, in the SF phase, the order 
parameter takes the form 



|Ao| = y/Eh{2n + Eh), 

implying a thermodynamic potential, ^ —i'{p + 
£'b/2)^, independent of the chemical potential difference, 
h. As a result, one may infer that the magnetization, 
m = — dfi/d/i, is zero. For the uniform condensate, 
the SF phase always involves a balanced population of 
fermions. Drawing on these results, let us now comment 
on the implications for the phase diagram of the system 
for, respectively, fixed chemical potential and fixed par- 
ticle number. 



A. Fixed chemical potentials 
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FIG. 1: The phase diagram shown as a function of 
{h/E^^, fj,/ E\j) for three different values of the mass ratio, r. 
The diagrams were constructed assuming a uniform order pa- 
rameter, neglecting the potential for inhomogeneous phase 
formation. The solid lines represent continuous phase bound- 
aries, while the dashed lines denote first order transitions into 
the balanced SF phase. 



When the chemical potentials, fi and h, are held 
constant, minimization of the thermodynamic potential 
leads to the phase diagram depicted in Fig. [TJ The equal 
mass case is consistent with the result of Tempore ct al. 
[23| . For II smaller than either the molecular binding en- 
ergy per particle, — £'b/2, or the chemical potential shift 
associated with the majority species, —h, (corresponding 
to an empty Fermi surface), the equilibrium phase hosts 
no particles (the "zero-particle" state, ZP). On increasing 
the chemical potential, fj,, a second order phase transition 
into either a balanced SF, or a fully-polarized (FP) nor- 
mal phase occurs. The transfers from the zero particle 



phase to a FP normal phase, and from a FP phase to 
a partially-polarized (PP) normal phase are both con- 
tinuous. The phases have boundaries where the Fermi 
surface shrinks to zero at /i = — /i and fi = +h respec- 
tively (for h > 0). At fixed i?b, an increase in chemical 
potential, fj,, leads to an increase in the order parameter 
of the balanced superfluid system, |Ao| oc ^y^, and an 
attendant increase in the critical h required to destroy 
the condensate. The phase transition from the normal 
state, both FP and PP, into the SF is first order. 

As the ratio of masses is increased, as shown in Fig. [T] 
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TABLE I: Summary detailing the loci of phase boundaries for ^/Eh as a function of h/Eh and R = [r — 1) / [r -\- 1). Results 
labeled {-h) are found in, and are relevant for, Sec. IIVI 

pp.ppLO* 1 + {h/E^)R ± yi + 2{h/E^)R -R^- 2{h/E^)R^ 



i?2 



FFLO-SF* and PP-SF ^ ^ E.)R - R^ ^ ^ - R^)^l + ^h/ E,)R) 



FP-SF 



2K^ 
V2{h/E^) + 



V2 - 2V1 - R 
ZP-FP ±{h/Eb) 

FP-FFLO* and FP-PP T{h/Ey,) 

ZP-SF -1/2 




FIG. 2: The phase diagram as a function of m/n and inter- 
action strength, Eh/Ep, of the two-dimensional system with 
fixed majority and minority particle densities for three dif- 
ferent mass ratios, r. The SF phase (dotted line) is the line 
of zero population imbalance. Between the balanced SF and 
PP/FP phase (dashed line) lies a region of phase separation 
(PS). 



B. Fixed number densities 



In the canonical ensemble, where the number densities 
n and m are held constant, the chemical potentials, fi 
and /i, must be inferred self-consistently. In this case, a 
first order transition in the {fj,/ E^, h/ E^) phase diagram 
(Fig. [1]) implies phase separation (PS) Q in the (n, m) 
phase diagram. At each point along the PP-SF phase 
boundary in {fx/ E\,,h/ Eb) one can evaluate the corre- 
sponding polarization and total number density. From 
this result, one can infer the boundaries between the nor- 
mal and phase separated regions as functions of E^/E-p 
and polarization, m/n. Here we have defined a "Fermi 
energy" scale Ep = n/v, where ly = m-R,/'2,n denotes the 
constant two-dimensional density of states of the reduced 
mass system. The resulting phase diagram is shown in 

Fig.m 

As expected, in the BEC limit of large E^/Ep, one 
finds phase separation, with the development of a con- 
densate of tightly-bound molecular pairs coexisting with 
a FP phase containing excess fermions. The phase di- 
agram shows that this behavior persists into the weak 
couphng BCS hmit, with the system phase separating 
into a balanced SF phase (i.e. with m/n = 0), and the 
excess particles forming a non-interacting PP Fermi gas. 
In the BCS limit of weak pairing, a small population im- 
balance is sufficient to destroy pairing and enter the PP 
normal phase region. 



on the side r > 1, the phase diagram becomes skewed. 
This can be understood by tracking the locus of the line 
where the Fermi surfaces of the two species are perfectly 
matched, approximately along the center of the SF phase. 
The central superfluid locus is ii/h = 1/R, which is con- 
sistent with the skew. Superfluidity is therefore more 
favorable if the "light" species has a greater chemical po- 
tential than the "heavy" species. 



When the species have unequal masses, the phase di- 
agram is skewed, similar to the fixed chemical poten- 
tial case in Sec. IIII Al If there is a mass imbalance 
then the Fermi energy of each spin species scales as 
E-p_cr = n^/va = T^in + crm)/m„ cx \/m„ implying that 
it is energetically more beneficial for the "heavy" rather 
than "light" particles to be in the normal state. There- 
fore, at a given mass imbalance, the phase diagram loses 
its symmetry in m/n and superfluidity is favored if the 
"lighter" species is in excess whereas the normal state is 
favored if the "heavy" species is in excess. 
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IV. INHOMOGENEOUS SUPERFLUID 

With the properties of the uniform SF phase in 
place, we now turn to the question of inhomogeneous 
phase formation. To characterize the nature of the PP- 
FFLO transition, we adopt two methodologies: firstly, in 
Sec. IIV A[ we will develop a Ginzburg-Landau expansion 
of the action to explore the locus of putative continuous 
transitions from the normal PP phase into the inhomo- 
geneous FFLO phase. Secondly, in Sec. IIVBI we will 
assess the validity of the Ginzburg-Landau expansion by 
investigating the global minimum of the thermodynamic 
potential for a mean-field order parameter field involving 
a single wave vector. Using these results, we will infer the 
phase diagram of a system with fixed chemical potentials 
in Sec. IIV C[ and fixed particle densities in Sec. IIVDI 



A. Ginzburg-Landau theory 

With the Ansatz that the transition from the normal 
to condensed phase is continuous, close to the transition 
we may expand the action in fluctuations, |Aq|. In doing 
so, one obtains 



•S'cff = y^QqlAql 



0(|A|4) 



(7) 



where 



E 



2ek + Eh ^k-q/2,T + 



denotes the static pair susceptibility. The locus of con- 
tinuous transitions may be determined from the value 
of q at which is both minimized and passes through 
zero. Within the condensed phase, higher order terms 
in Aq determine the crystalline structure of the FFLO 
state 

The corresponding phase boundary then translates to 
the largest allowable chemical potential shift, h, which 
occurs when the Fermi surfaces just touch but do not 
cross [HH. From this condition, one finds a phase bound- 
ary along the line, 



Eb \Eh 



(8) 



Minimizing aq with respect to |q|, one obtains the further 
condition eq = q'^/{2m^) = 2Eh/{l — R^). Measured 
in units of the Fermi momentum of the reduced mass 
system, this translates to a wave vector. 



M 

k-p 



I Eh {mi +mi) 
Ef 2mR 



(9) 



where fcp 



2mY{EF and, inverting Eqn. ([8]), Eh 



r = 2 



6 fi/Eh 




0.2 

0.15 

\^Q\/Eh 0.1 
0.05 




6 fi/Eh 



FIG. 3: Variation of the order parameter | Aql/i^b with chem- 
ical potential difference, h/Eh, and chemical potential, fi/Eh- 
The upper panel is at a mass ratio r = 2, and the lower panel 
at equal masses, r = 1. The thin black lines trace out the 
|AQ|/_Eb variation, found by minimizing the thermodynamic 
potential Eqn. for several different chemical potentials. 
The thick black line is the locus of the second order transi- 
tion predicted by Ginzburg-Landau theory. 



limit. Eh <C Ep, so that, at equal masses. Eh = h? /2ii 
giving |q| = 2h/vF, where vy is the Fermi velocity, 
agreeing with the findings of Burkhardt and Rainer 
[63 |. Shimahara [tO], and Combescot and Mora 52 1. In 
the same limit, the pair susceptibility takes the form, 
5R(ln(l + \Jl — (|q|tiF/2/i)^)), collapsing to that fomid in 
previous works. 



hR ± — /i^)(l — R^)- In the weak coupling 



B. FFLO instability phase boundaries 

To assess whether the transition from the PP phase 
to the FFLO phase is really continuous, one can instead 
minimize the thermodynamic potential Eqn. ([1]) with re- 
spect to the wave vector Q and the mean-field value of the 
order parameter Aq. For several values of chemical po- 
tential, and two different mass ratios r = 1 and r ~ 2, 
numerical minimization of the thermodynamic potential 
confirms that the order parameter changes continuously 
(see Fig.|3|), falling to zero along a line of instability. The 
locus of the transition also agrees with that obtained from 
the Ginzburg-Landau expansion in Sec. IIV A1 This result 
is in accord with that found in Ref. [l^l in the weak cou- 
pling limit of the equal mass system, and shows that the 
transition remains continuous across the entire range of 
the FFLO phase. 

We are now in a position to evaluate all phase bound- 
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aries associated with the FFLO instabiUty. The agree- 
ment described above between Ginzburg-Landau theory 
and direct minimization allows us to use the analytic 
Ginzburg-Landau boundary between the PP and FFLO 
phases. The minimum in the thermodynamic potential 
that gives rise to the FFLO phase is shallow relative to 
that of the SF phase. We are therefore able to approxi- 
mate the actual FFLO-SF phase boundary by the Q = 
result for the PP-SF boundary described in Sec. IIIII A 
summary of the phase boundaries is shown in TablelH the 
additional boundaries due to the presence of the FFLO 
phase are labeled (if). As the extent of the SF region is 
only reduced by the presence of the FFLO phase, the SF 
is balanced, as was shown for the Q = study in Sec. IIIII 



C. Fixed chemical potentials 

Let us now apply these results to the problem of a 
uniform system with constant chemical potentials. The 
corresponding phase diagram is shown in Fig. [D While 
the general topology of the phase diagram mirrors that 
discussed in Sec. IIII Al the transition to the balanced SF 
phase is preempted by the formation of an inhomoge- 
neous FFLO phase. The FFLO instability occurs mainly 
on the PP side of the PP-SF phase boundary of the uni- 
form condensate shown in Fig. [T] with just a small in- 
trusion on the balanced SF side. The FFLO instability 
does not occur within the FP state as there are no mi- 
nority state particles with which to pair. The FFLO-PP 
boundary is second order, while the FFLO-SF boundary 
is first order. 

The FFLO-PP phase boundary terminates at the SF 
phase for small mass ratios and at the FP phase for large 
mass ratios on the side of the majority "heavy" species. 
The movement of the boundary terminus with increasing 
mass ratio r is in the opposite direction on the majority 
"light" species side - it moves further up the SF phase 
boundary. The special mass ratio where it terminates 
at the SF-FP phase boundary on the majority "heavy" 
species side is at Tc = (1 -I- V5)/2. 

The thermodynamic potential variation is also shown 
in Fig. 3] at four different points (a, b, c, d) for r = 1. 
Since the wave vector dependence of the thermodynamic 
potential enters through the order parameter, in both 
the PP (a) and FP (d) normal phases the minimum is 
|Q|-indcpendent. At the highlighted FFLO phase point 
(b), the global minimum lies at |Aq| w 0.2i?b with 
IQI ~ 2vmR^Eb, while a local minimum also develops 
at |Aq| w 3.8£'b with |Q| = corresponding to the pu- 
tative uniform SF phase. At the highlighted SF point (c), 
the global minimum lies at |Aq| « 3.8E\^ and |Q| = 0. 

D. Fixed number densities 

Let us now address the implications of the phase dia- 
gram for a spatially uniform system held at fixed num- 



ber densities. Obtaining the corresponding density, n, 
and magnetization, m, gives the phase diagram shown in 
Fig. [5] Once again, the topology of the phase diagram 
mirrors that discussed for the homogeneous condensate 
in Sec. IIII Bl However, between the phase separated SF 
phase and normal phase, the system exhibits an inho- 
mogeneous FFLO phase over a wide region of the phase 
diagram. 

In the weak coupling BCS limit, even a small popu- 
lation imbalance is sufficient to enter the FFLO phase 
region. We note that in a population balanced system 
the Fermi momenta of the populations are equal so no 
shift of the Fermi surfaces is required to form Goopcr 
pairs and a modulated phase is not seen. 

The effects of the moving PP-FFLO phase boundary 
terminus, described in Sec. IIV Gl are also apparent. For 
equal masses, the FFLO phase never meets the FP nor- 
mal state. For mass ratios in excess of rc = (1 + \/5)/2, 
the FFLO phase meets the FP state on the majority 
"heavy" species side, but is further from the FP state 
on the "light" species side. For high mass ratios, this 
is evidenced by the much broader FFLO region on the 
heavy species side. 

To conclude this section, it is interesting to compare 
the phase diagram of the ultracold atom system with con- 
tact interaction and the problem of electron-hole bilayers 
with long-ranged Goulomb interaction. In particular, we 
focus our discussion on the study in Rcf. [13 of GaAs 
bilayers where the mass ratio r ~ 4.3. In this case, it is 
more natural to characterize the strength of interaction 
by = ro/op, where rg = denotes the interparti- 

cle spacing, and uq is the effective Bohr radius of the two- 
body bound state. The latter is related to the dimension- 
less ratio E'a/E-p through the relation, E\^/Ey = O.BSlrg . 
As a result, we find that the system enters the BGS phase 
with the appearance of FFLO phase behavior for 7's val- 
ues of ca. 1.5 (4) compared with that found for the un- 
screened electron-hole bilayer of Tg ^ 1.5 (16) for the 
"light" ("heavy") species. More qualitatively, in both 
cases, the systems show a preference towards the super- 
fluid phase when the "light" species is in excess, and 
the normal phase when the "heavy" species is in excess. 
Although the topology of the phase diagram is quantita- 
tively the same, two significant differences appear. The 
first is that, with the electron-hole bilayer, the FFLO-SF 
phase boundary on the "heavy" species side penetrates 
further into the BEG regime than in the ultracold atomic 
gas. The second difference is that, with the electron-hole 
bilayer, the FFLO region existed from the normal phase 
to TO = 0, and no phase separation between FFLO and SF 
was seen, except for the deep in the BEG regime. How- 
ever in the ultracold atomic gas, phase separation of the 
SF was seen into a balanced SF and a FFLO phase. Both 
of these differences indicate that, with the electron-hole 
bilayer, the FFLO phase was more stable relative to the 
SF than in the ultracold atomic gas. This could be due to 
the long-range forces that act in the electron-hole bilayer 
whereas the ultracold atomic gas experiences only con- 
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FIG. 4: The left-hand column shows the phase diagram in {^./ Eh,h/ Eh) at three different mass ratios r. The solid lines 
represent second order phase boundaries, the dashed line denotes first order phase transitions into the balanced SF phase. The 
FFLO phase is highlighted in grey. The trajectories followed in the sample traps in Sec. |V] are shown by thin dashed lines. 
The central column of diagrams focus more closely on the topology of the phase diagram close to the tricritical point region for 
h/Eh < 0. Thermodynamic potential surfaces are shown in the right-hand column; the darker the more negative and so more 
favorable; label (a) is a PP normal state, (b) a FFLO state, (c) a SF state, and (d) a FP state. Different plots have different 
shading calibrations, (unequally spaced) contour lines are also shown. 



tact forces that would favor formation of tightly-bound 
BEC pairs. 



V. HARMONICALLY TRAPPED SYSTEM 

Finally, focusing on applications to ultracold atomic 
gases, we now address the influence of the trap geometry 
on the phase behavior. Here we make use of the local 
density approximation in which the chemical potential 
of both species, ^^(R) = /i^ — V^(R), are renormalized 
by the same local trapping potential 1^(R), the chemi- 
cal potential difference, h, remains fixed across the trap. 
Moreover, we further assume that the spatial coordinates 
are rcscalcd to ensure a spherically symmetric trapping 
potential, ^(R) = luR^2. Although there IS some ex- 
perimental evidence [tH . [T^I that the local density ap- 
proximation might not be valid p^ . [t^ in some setups, 
we believe that its application here will correctly address 
the qualitative phase structure. 

To identify the phases present, one may consider a tra- 
jectory of changing fi with constant h and r in the phase 
diagram of fixed chemical potentials. To find the to- 
tal magnetization and number of particles in the trap, 
one may make use of the local relations m = —dQ/dh 
(Eqn. (12)) and n = — dfi/d/i (Eqn. ([3])) respectively, and 
then integrate over the trap. All trajectories will end up. 



at large enough radius, in the ZP regime, which is the 
edge of the particle distribution. 

The profiles in four sample traps arc shown in Fig. [51 
which follow trajectories highlighted in Fig. 21 The first 
three have species with equal masses, r = 1. At zero pop- 
ulation imbalance only the SF state is observed. With a 
population imbalance, firstly there is a central balanced 
SF region surrounded by a ring of FP majority spin par- 
ticles. On increasing the population imbalance yet fur- 
ther, between the ring of FP particles and the central 
SF, an FFLO instability adjacent to a PP region is seen. 
The first order transition between the SF and FFLO re- 
gion (and FP state) leads to a discontinuity in density 
and polarization. The second order transitions between 
FFLO, PP and FP states have continuous variation of 
density and polarizability but discontinuous changes in 
their gradients. 

The final profile in Fig.|6lis at an unequal mass ratio, 
r = 4. The inclusion of mass imbalance causes the SF re- 
gion in Fig.[4lto be biased towards the "lighter" species. 
This means that it is possible to have a ring of superflu- 
idity remote from the trap center, or an isolated ring of 
FFLO instability not at the center and no accompanying 
SF region. When there are two rings of normal phase 
bounding the SF they may either both be the "heavy" 
particle normal phase if we are crossing the extrusion of 
the FFLO phase, or alternatively one might be "light and 
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r 




FIG. 5: The phase diagram in {m/n, Eb/Ep) of the two-dimensional system with fixed majority and minority particle densities 
for three different mass ratios r. The SF phase (dotted line) is represented by the line of zero population imbalance m/n — 0. 
The PP phase is separated from the FFLO phase (highlighted in grey) by a second order phase boundary (solid line). In the 
phase separated (PS) region bordered by the dashed line, the system separates into a balanced superfluid (SF) and either a 
normal phase or, depending on the composition, an inhomogeneous FFLO phase. The right-hand column of graphs focus on 
the terminus of the PP-FFLO phase boundary on the side of negative polarization for the three featured mass ratios. 




R^f^Wh R^^/E^ R^uj/Eb i?vW^ 



FIG. 6: (Color online) Radial density profiles of four systems in identical harmonic traps with different values of population 
imbalance, h, and mass ratio, r. The solid line shows the radial density based on the primary y-sxis. The dashed line shows 
the local population imbalance based on the secondary j/-axis. The band shows where different phases exist with colors labeled 
on the right-hand side of the figure. With changing radial coordinate, the separate panels are associated with the trajectories 
highlighted in the {h/ E^, fi/ E^,) phase diagram shown in Fig. |4] 



the other "heavy" if traversing right across the skewed SF 
phase. In the latter case, shown in Fig. [6l the species fa- 
vored by the chemical potential shift dominates at the 
outside of the trap. At the center of the trap, the normal 
state is of the "heavy" species as superfluidity favors the 
"lighter" species. 



VI. CONCLUSIONS 



We have derived an analytic expression for the thermo- 
dynamic potential of a two-dimensional two-component 
atomic Fermi gas in the mean-field approximation with 
population imbalance and general mass ratio at zero tem- 
perature. A complementary Ginzburg-Landau analysis 
was used to examine the PP-FFLO transition. Ana- 
lytical expressions for the phase boundaries separating 
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normal and inhomogeneous superfluid phases have been 
obtained while the properties of the FFLO phase have 
been addressed numerically. Within the mean-field ap- 
proximation, the SF phase does not sustain a population 
imbalance. The region of FFLO instability exhibits a 
second order phase transition from the PP normal phase, 
and first order phase transition into the balanced SF. In 
the BCS limit, a small population imbalance is sufficient 
to destroy pairing. In the BEG limit, there is phase sepa- 
ration between tightly-bound molecules and a FP normal 
phase. If there is a mass imbalance, the SF phase is fa- 
vored if the majority particles are the "lighter" species, 
while the polarized normal state is favored if the "heavy" 
species are in excess. 

A trapped geometry leads to a rich range of possible 



density profiles. If there is no mass imbalance, a SF phase 
is seen at the trap center surrounded by a PP followed 
by a FP normal phase of the majority spin species. If 
there is mass imbalance, then a ring of the SF and/or 
the FFLO state could be seen bordered both inside and 
outside by either species of normal phase particles. 
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